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ON STABILITY OF FREE LAMINAR BOUNDARY 
LAYER BETWEEN PARALLEL STREAMS 

By Martin Lessen 


SUMMARY 

An analysis and calculations on the stability of the free 
laminar boundary layer between parallel streams were made 
for an incompressible fluid using the Tollmien-Schlichting 
theory of small disturbances. Because the boundary condi- 
tions are at infinity, two solutions of the Orr-Sommerfeld 
stability equation need not be considered, and the remaining 
two solutions are exponential in character at the infinite 
boundaries. 

The solution of the stability equation is obtained in powers 
of (— ifaR ), where a is disturbance wave number and R is 
Reynolds number. With an asymptotic solution as a start, 
ike stability equation is numerically integrated. The eigen- 
value problem between a, R, and the disturbance phase velocity c 
was thus explored by trial-and-error process. 

The calculations show that the flow is unstable except for 
very low Reynolds numbers. The regions of stability and 
instability in the a,R-plane were checked by obtaining damped 
and amplified solutions on opposite sides of the neutrally 
stable solution. 

INTRODUCTION 

Some of the classical problems that have been of interest 
to many investigators in the field of the stability of parallel 
flows are the stability of Couette type motion, Poiseuille type 
motion, and boundary-layer flows. Couette flow between 
rotating cylinders was successfully treated by G. I. Taylor 
(reference 1). The problem of plane Couette flow, however, 
has not been decisively settled. Plane Poiseuille flow was 
treated by Heisenberg (reference 2), Pekeris (references 3 
and 4), Goldstein, and Lin (references 5 and 6) and highly 
controversial results were obtained. Poiseuille motion in a 
circular pipe was treated by Sexl (references 7 and 8) who, 
after an incomplete investigation, concluded that the flow 
was stable. The problem of boundary-layer stability was 
treated by Tollmien (references 9 to 11), Schlichting (refer- 
ences 12 to 15), Lin, and others, and an experimental verifi- 
cation of the results was obtained by Schubauer and Skram- 
stad (reference 16). 

In all the previous investigations, the problem of the sta- 
bility of parallel flows was treated for a case in which there 
was at least one solid boundary. In the investigation dis- 
cussed herein, a flow that has no solid boundaries is treated, 
the case of the stability of the laminar, free boundary layer 
between two parallel streams of semi-infinite extent in plane 
flow. In the analysis, the fluid medium is considered incom- 


pressible. After the steady-state flow configuration is de- 
termined, the stability of that flow configuration is investi- 
gated. The solution is carried out for only the case of flow 
in which one of the streams is considered at rest. 

Most of the analysis for this report was done before June 
1948 in the form of a doctorate thesis in the Mechanical 
Engineering Department at Massachusetts Institute of 
Technology. Acknowledgment is made to Prof. W. R. Haw- 
thorne of the Mechanical Engineering Department at 
M. I. T. for his general supervision of the work, to Prof. C. G. 
Lin of the Mathematics Department at M. I. T. for his 
valuable assistance and suggestions, to Dr. L. H. Thomas 
of the Watson Scientific Computing Laboratory at Columbia 
University for his assistance in setting up the problem for 
machine solution, and to the International Business Machines 
Corp. for donating the use of the Selective Sequence Elec- 
tronic Calculator for the numerous calculations. 

SYMBOLS 

The following symbols are used in the theoretical develop- 
ment: 

x positional coordinate in direction of principal flow 

y positional coordinate perpendicular to principal 

flow 

t time coordinate 

u velocity component in ^-direction 

v velocity component in y-direction 

stream function 
f vorticity 

p mass density of fluid medium 

n absolute viscosity of fluid medium 

v= nip kinematic viscosity of fluid medium 

S characteristic measure of boundary-layer thickness 

Ul free-stream velocity of one of parallel streams 

(taken as characteristic velocity) 

U 2 free-stream velocity of other parallel stream, 

U{>U 2 

R—sUi/y Reynolds number of boundary-layer flow 
i]= v! fvxjUi 

f(jj) function defining form of boundary-layer stream 

function for time-independent flow 
x*=xfS dimensionless positional coordinate in direction of 
principal flow 

y*=yfS dim ensionless positional coordinate perpendicular 

to principal flow 
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dimensionless time coordinate 


u*—u!Ui dimensionless velocity component in redirection 
v*—v/Ui dimensionless velocity component in y*-direction 
f* dimensionless stream function 

r^Hy*)+4>(y*)e ia(x, ~ cl,) 

${y*) steady-state part of dimensionless stream function 

amplitude of disturbance part of dimensionless 
stream function 


* V )= 


/' 0 ?)= 


a 

c 


c, 


Erjc z — 


h\, ki, ks 

Im 

exp 


dcj> 

df* 

df 

di) 

wave number of disturbance (always real) 
eigenvalue of a when R— * °° 
dimensionless phase velocity of disturbance 
eigenvalue of c when B-* °° 

V=T 

j* dt 

constants of integration 
imaginary part of 

base of Napierian logarithmic system e raised to 
power in parentheses following exp 


and equation (1) simplifies to 

//"+2/'"=0 


( 3 ) 



From the physical flow configuration (fig. 1), the boundary 
conditions are 

y— S-+CO u-*Ui 

y-^—co u-*U t 


THEORETICAL CONSIDERATIONS 


or, in terms of the variables of equation (3), 


STEADY-STATE. LAMINAR, BOUNDARY-LAYER FLOW BETWEEN PARALLEL 

STREAMS 

The equations for the steady-state incompressible boun- 
dary-layer flow with no body forces and no pressure gradient 
over the flow field are 


bu , bu d 2 w 

U bi +V bi =v w 

(1) 

bu , br 
bx + by~ 

(2) 


From continuity considerations (equation (2)), a stream 
function f can be so introduced that 


77 — >-}- CO 

17— > — 00 

Because: equation (3) is of the third order, a third boundary 
condition must be. specified. The third boundary condition 
is arbitrarily selected as 

V = Va J= 0 

It is of interest to discuss the asymptotic behavior of the 
boundary-layer equation (equation (3)). When 

ij— >00 


1 

J U: 



bf 


bx 


If the form of ^ be specified as 


/'(*)-> 1 

and 

Therefore, as a first approximation for largo positive values 
of rj, equation (3) can be written as 

Vf'+2f"^0 


1 P= fx U- f(ri) 


and it therefore follows that 


then 


V=vl^vx/U 1 


u=UJ'( n ) 


l-f^kiErjc(£) (4) 

Similar!}', for large negative values of n, 



JvUijxXnf-f) 



(5) 
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Equation (5) degenerates for U 2 jUi= 0. For tliat case, if 
it is assumed that when 




CO 


Equation (8) therefore reduces to 
/(»)[!■ -^ = 0 


o 

and 

a 

then 

fto-+k* v 

It is also of interest to discuss the position of the boundary 
between the fluid of both streams. This boundary must be 
a streamline and also pass through the point (x=0, y=0). 
From the form of the stream function yp, it can then be seen 
that ip—Q is the boundary streamline. In order for this 
condition to hold at any value of x, it is necessary that/(i?) =0 
along the streamline. The boundary streamline is therefore 


The significance of the v» hoe can also be demonstrated 
from momentum considerations. If, for a fixed value of x, 
the boundary streamline is located at y t , the momentum 
principle can be formulated as follows: 

pu(Ui— v^dy+j" 1 pu{U t —u)dy—Q (6) 

If/ and ij are then substituted into equation (6) , the following 
equation is obtained: 



From the asymptotic behavior of / (equations (4) and (5)) 

/(«)[!— /'(») 1=0 

l 

/(-«)[§-/'(- «)]=° 

Therefore, 

/(’ii)=0 

and jji=i 7 o, as previously demonstrated. 


STABILITY OF LAMINAE BOUNDARY LAYEE BETWEEN PARALLEL STREAMS 

The criterion used in the investigation of the stability of 
the laminar -flow configuration is the behavior with respect 
to time of a small periodic disturbance introduced into the 
flow field. If the disturbance is damped with respect to 
time, the flow configuration is said to be stable with respect 
to that disturbance. If the disturbance grows with respect 
to time, the flow' is said to be unstable. If the disturbance 
remains unchanged with time, the flow is considered as 
neutrally stable. 

In reference 17, Squire demonstrates that two-dimensional 
disturbances are more destabilizing than three-dimensional 
disturbances for two-dimensional parallel flows. Only two- 
dimensional disturbances are therefore considered in this 
discussion. 

The equation of vortieity for an incompressible viscous 
fluid in plane flow can be stated as 


JJ /'(i -f)dv +/^ f (vr f ) dv = 0 

where ij t is the boundary between the two streams. Then 
/(«)-(! -^) Jiv 0 -§/(--) - (W V = 0 (7) 


where 


However, 


/(/')*«? iJ 


I'Af 


dr du 
* dz dy 





Af-a+S 

s dr 2 dy s 


( 9 ) 


From equation (3) 


Therefore, 


//"=- 2/' 


f(jyd v =ff'+2j" 

When the preceding relation is substituted into equation (7), 

/(”)—(_! — £r) /0 7 0— 

2 /"(»)+ 2/"(-»)=0 ( 8 ) 

However, 

/"(-)=/"(- <*0 = 0 


Now, let f , u, and u consist of steady-state (time-independent) 
and disturbance (time-dependent) parts. Then 

r=?(z,y)+r'(z>y,f) 

u—u(x,y)+u'(x,y,t) 

v=~c(x,y)+v'(x,y.t) 


where 


r'«r ti'«u v'«d 


If the expression for f, u, and v are substituted in equation 
(9), the steady-state terms cancel out as satisfying them- 
selves. Because the disturbance is considered very small 
compared with the steady-state flow, all terms nonlinear in 
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Equation (13) then becomes 


the disturbance are neglected. Equation (9) then becomes 

It is assumed that u' and v' are of the same order of 
magnitude and that and are of the same order of 
magnitude. If it is now considered that the steady-state 
flow is of the boundary-layer variety ’and that the principal 
direction of flow is the r-direction, then 

v«u and 7x«T* 

Equation (9) then reduces to 

r 7 « + «£ '*+ v'$y— vAf ' 

If a stream function of the form 

/*)+0(yV aC **- e “ > 


( 10 ) 


( 11 ) 


is introduced into equation (10) after the dimensionless 
variables u* t v*, x*, and y* are substituted for the corre- 
sponding dimensional variables, the Orr-Sommerf eld equation 
is then obtained (references 18 and 19). 

($'-c)(<J>" - a 2 <t>)-&"<t>= (fa v -2 a 2 <f>"+ a*4>) (1 2) 

where 

T f ^ ft! V 

* = ^ ?==w (v) 

A more rigorous and detailed derivation of the Orr-Sommer- 
feld equation may be found in reference 20. 


y —n 


If the characteristic length 8 is set equal to y -jj> then 
and 


dy* 3 dy 3 9 

The Orr-Sommerfeld equation for the laminar boundary 
layer therefore becomes 


</-c)(<r-«V)-/"'*=-^(* ,T - 2aV'+«V) (13) 

For boundary conditions on equation (13), it suffices to 
specify that <j> is bounded throughout the flow field. 

Because equation (13) is of the fourth order, a set of four 
linearly independent solutions to the equation exists. 

<t> = C 1 fa + C ifa + C&fa + C t fa 

The nature of the solutions fa and fa can be investigated by 
introducing the following transformation into equation (13): 

<t>=exp{f gdy) 


~ [g t +69V+Zg' t +^gg"+ 

g'"-2 cdigt+g')-**] 

If g is expanded as follows: 

j_ i 

g=(aK) 2 g 0 -{-gi-h(aR) 2 . . . 

the following set of equations is then obtained by equating 
like powers of aR: 

(f'—c) ga=—igo A 

(f - C) (go + 2mi) = —i (tyoOi + 6yoV) 


The successive approximations can then be obtained with- 
out integration 

Po= ± jiif—c) 


9 1 2 0 O 


Finally, 


fa =ss(f— c) * exp f -\/i <xR(f—c) dyj 

faez(j'-c)" exp [+ j Vi ceR(f—c) dyj 
It should be noted that when 

TJ — > OO (j)g — > CO 

1j — > —|— OO ^ — > 00 

Because fa and fa are each unbounded somewhere in the 
flow field, they cannot be considered in the solution of equa- 
tion (13). 

The quantities fa and fa can be solved for in the following 
manner: 

Let 


*=S(-a?)> 


(14) 


If equation (14) is substituted into equation (13) and like 


orders of I 
tions is obtained: 


( <Zr) 


are equated, the following system of cqua- 


( 15 ) 

For k=0, equation (15), after being rearranged, becomes 

*<o>"-(W^) ,»<<» = 0 (16) 
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The asymptotic behavior of 0® 
Where 


ij — »± <» 


can now be investigated. 


Therefore, 


_f 1 

f~c 


*0 


Because the solution must remain bounded, 

<?i=0 when »+“ 
t?i=0 r \— * — oo 

In the region where <f>® ) is exponential, <f> a) is also exponential, 
as can be shown by substitution into equation (15). It can 
thus be demonstrated that $ a) is also exponential in the 
aforementioned region. The proper boundary conditions 
on <(> are therefore 


ij— > — oo <j>'—^acj) 

1J— >-f-00 cj)'— > — a4> 


From the boundary conditions and the general form of the 
solution, the secular relation formally stating the eigenvalue 
problem can be easily obtained. 


F(a,c,R) 


0i'(— 00 )— <*&(— 00 ) 

<t> /(+ 00 )+ °°) 


°°) — a0j(— ~) 
<£/(+ °°) + «02(+ 03 ) 


= 0 
(17) 


From equation (17), it can be seen that only particular com- 
binations of the parameters a, c, and R allow the boundary 
conditions on 4> to be satisfied. 

For a rigorous treatment of the asymptotic expansions for 
02, 03, and 04, see reference 21. 


SOLUTION OF BOUNDARY LAYER AND ORR-SOMMERFELD 
EQUATIONS FOR LAMINAR BOUNDARY LAYER BETWEEN 
PARALLEL STREAMS 

The general plan of solution of the boundary layer and 
Orr-Sommerfeld equations is to cany the analytical methods 
only far enough to reduce the problem to a routine numerical 
solution. In the problem under consideration, it is neces- 
sary to investigate the asymptotic behavior of the solutions 
so that boundary conditions at finite values of the independ- 
ent variable are known in terms of the infinite boundary 
conditions. After boundary conditions at finite points are 
obtained, the solution can be continued by ordinary numeri- 
cal methods. 

ASYMPTOTIC BEHAVIOR OP STEADY-STATE BOUNDARY-LAYER EQUATION 

The following transformation of equation (3),/f"+2/'"=0, 
is permissible without altering the form of the equation: 

s(ar)=^/(’7) 

x=at}-{-b 

956646 — 51 3S 


If the discussion is limited to the case where U 2 /Ui=0, it is 
then possible to expand a solution for the transformed equa- 
tion as follows: 

As A r 

2 = .Ao - ! - A.],e 2 -f- A-ie x -r A$e 2 + . . . 

The foregoing expansion is convergent for negative values 
of x. If a numerical solution is started from inside the area 
of convergence of the series solution, the asymptotic value of 
g'(co) can be evaluated. Because the asymptotic value of 
/'( co ) is ij the scaling factor a is given by 

1 

a= — 

vV( ro ) 

The scaling factor b can be fixed by designating 

/W= 0 

The asymptotic form of/ is then given by 

-an 

/=a 0 +a ie 2 +a.2e a, i+a 3 e 2 + . . . 
and, because/' (starting) </, 


-Ml -ar, 

=M 2 n + b 2 e‘*+b 3 e* 

J — c 

The various coefficients of the expansion can easily be 
evaluated by elementary methods. 

ASYMPTOTIC BEHAVIOR OF EXPANDED ORR-SOMMERFELD EQUATION 

If the asymptotic form of/'"/ (/' — c) is inserted into equa- 
tion (16), the following equation is obtained: 

( I 3 „ \ 

+ b 2 e°i+b s e 2 +...,) 0«> = 0 
and the form of 0 CO) follows: 

0® )=e c^ +h ^ a ^ a > +h2ae (a + ^ + h 3 J a+Ta ' ) ’’+ . . . (18) 


If <j> m is substituted into equation (15) to obtain the next- 
order term 0 (1> , the following equation is obtained: 

<j> iiy'-( a z+b 1 e Ta ’'+b 2 e°'+b 3 e Ta '+ . . . ) 0 (1) 


=l l e^ ai ~ 2a ^ +lie< a+a) '+l 3 e(- a+T ' I ' > ' 

and 

(oe-i-i- a)? (a-r-f-a)* 

0 (l > =e a '+hne^ -\-k 2l e^<*'+h 3l e^ 2 J + .. . (19) 


From equations (18) and (19), the boundary conditions 
can be obtained for integrating the stability equations 
starting at a finite 17. 

A tabulation of the values for the coefficients of the first 
few terms of the expansions is given in the appendix. 
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INTEGRATION OF BOUNDARY-LAYER EQUATION 

The boundary-layer equation is 

//"+ 2/'"=0 (3) 

or 

Differentiation once gives 

/ lv — |c ff"+m 

and twice, 

f=~\ (jf lv + 2 ff"+f' 2 ) 

Expansion in a Taylor’s series gives 

Av+w)=M+wf(v)+~f'W+^-f'^^ 

f(v +v>) =/'<*) + wf'iv) +■ ~f"(v) ./ lv M + 

.... -- ' 

f'(r,+w)=r(r,)+wf"( »?)+y/ lv W+yV T (n)+ 

■ * 

where . 

Sf(v)=f(v)-f(v~w) 

The values of /"'Oj+w), y iv (i?+w), / T (r+w), and 5j v (rj + w) 
can now be algebraically computed, as indicated and the 
integration carried forth for the next interval. The process 
can be started by evaluating 5/ v from the asymptotic form 
of the solution. 

INTEGRATION OF EXPANDED ORR-SOMMERFELD EQUATION 

Consider the equations for the first two terms of the 
expansion (equation (15)) of <j>: 

(/' - c)(4 > C0) " - a 2 4> (0 >) -/'"«*> (0> = 0 

{f -c)(«f> (1) = (<t> m !v - 2 a 2 ^ " + aV 0) ) 

Expansion of <£ (0) in a Taylor’s series gives 

<t> w (v+u>)= 4> {0) (v)+w<t> m '(v)+y ^ "(v)+~<P (0) '"(v)+ 

£ <f> (0,lr (v)+ ■ • • 

y '(v)+Y 0 CO)IV (n)+ 

^ ^« *•(?)+ — ■ 


The values of <f> (0) "(j? + w), # c<s "0? + ^ craiv (i?hiu), and 
6 ^ (0) Iv (t;+w) can now be computed in a maimer similar to 
the case of the boundary-layer equation and the integra- 
tion of 0 tO) can be carried forth. As before, 5<f> <0) lv can bo 
evaluated in starting. 

Expansion of 4> U) in a Taylor’s series gives 
4> li Kv+w)=^Hr,)+w^'(v)+Y 4> ,l), \v)+^ <£ fl> '"(n)+ 

^ l, 1 v W+^T * Wv (v)+^ 0 <l,T, (»)+ ■ • • 


Now 8 i <j£±(ii)=8<l> (l> (,r}-{-w ) — 

Therefore, 

5 2 <f> ® (i?) = <j > (1) (ij -f- w) + (q — w ) — 2 <£ (l) (ij) 

=^ < ‘ , -W+^^ tl,,v (»))+^^ a,T ^)+ • • • 

and ....... 

a V , )"(^)= W V« )l v W - [ -^ 0 a,vt( n ) + . . . 

Therefore, 

5 2 [t m (ri)~ ^W"(,)]«W*# n >"(l,) 


The quantity <£ 0) can be integrated as indicated in a step- 
by-step manner. From the asymptotic form of 4> (U , 

5 can be evaluated in starting. 


When 0 (o) and <£ (l) have been integrated to a largo positive 
value of ??, the boundary condition at 7= + <» can be imposed 
to evaluate B\ 

*“*" + (5?)*“ 

(^’(-j- <»)-(- Qf^(-(- Co) = Q 


Therefore, 

»_i r <f> w '(+”)+^ (1) (+”) i 

«U < ®.'(+“)+«^(+®)J 


( 20 ) 


Although the foregoing metliod of integration was used in 
the solution of the problem, it was not the most efficient 
scheme. Originally, the problem was programmed for 
integration along a rectangular path through the complex 
plane; complex integration was necessary because of the 
singularity at the point where f=c in equation (15). It 
later became apparent that the rectangular path was dis- 
advantageous and that a better choice was a path starting 
in the third quadrant of the complex plane and traveling in 
a straight line to a point on the positive, real axis (fig. 2). 
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Figcre 2.— Path of numerical integration of boundary-layer and disturbance equations in 

complex 11 -plane. 


Most of the original program could be used for the straight- 
line path. In order to expedite the solution, the equations 
were therefore integrated along the new path in the manner 
previously indicated. The interval of integration w was 
taken as 0.2+0.05 i. 

A more advantageous scheme for the straight-line path of 
integration of equation (15) is included herewith for complete- 
ness. 

52 

METHOD OF SOLUTION OF EIGENVALUE PROBLEM 

If sets of eigenvalues of a, c, and R are to have any physical 
significance, R must be real; a is taken real and c may be 
complex. The process then to be followed is: 

1. For a fixed value of a and a set of values of c, integrate 
the equations and obtain the corresponding values of R. 
These values of R are usually complex. 

2. For the value of R that is real, the corresponding c is 
the desired eigenvalue. 

3. Repeat the process until the secular relation (equation 
(17)) is explored. 

RESULTS AND DISCUSSION 

The secular relation (equation (17)) for the case of the 
‘free” boundary layer was explored for damped, neutrally 
stable, and amplified disturbances. The stable disturbance 
corresponds to Im e<0, the neutrally stable disturbance 
corresponds to Im c= 0, and the unstable disturbance cor- 
responds to Im c> 0 (equation (11)). The disturbance 
equations were solved for values of Im c= 0, Im c=0.05, 
Im c=0.10, and Im c—— 0.05 (table I). The results were 
then interpolated for Im R = 0 to give the neutral stability 



curve and curves of equivalent degrees of damping and 
amplification (table II and fig. 3). 

The curves obtained were continuous, which indicates a 
low effect of rounding errors in the numerical techniques 
employed. Some integrations were redone using an interval 
of half that used throughout the problem. The results of 
both integrations were the same to the fifth significant 
figure and indicated the low truncation error. Because the 
parameter of expansion of the eigenfunction (equation (14)) 
was (-11 aR) and only the first two terms of that expansion 
were used, the determination of the eigenvalues of a, c, and 
R was inaccurate at low values of aR. For the same reason. 



Reynolds number-, R 

Figure 3.— Curves o £ equivalent amplification and damping for free boundary layer. 


it was impossible to obtain the lower branch of the curve of 
neutral stability. For purposes of comparison with the 
free boundary layer, Lin's curve of neutral stability for 
Blasius type flow is included in figure 3, (See references 5 
and 6.) 

A comparison of the stability characteristics of the free 
boundary layer with Blasius type flow reveals certain basic 
differences. First, the free boundary layer is unstable for 
some disturbances at an infinite Reynolds number (tanta- 
mount to "in viscid” flow), whereas the Blasius type flow 
is completely stable in that range. This type of instability 
is caused by the presence of a point of inflection in the 
velocity profile. Some implications of the presence of a 
point of inflection in the velocity profile were first investigated 
■ by Rayleigh (reference 22). Another salient difference 
between the two types of flow is the presence of a highly 
oscillatory solution (<&) in the general solution for the dis- 
turbance in Blasius type flow, whereas no such solution 
exists for the free boundary layer. This difference arises 
from the doubly infinite boundary conditions for the free 
boundary layer as opposed to the singly infinite boundary 
conditions for Blasius flow. The foregoing is apparent from 
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the asymptotic forms of the solutions fa and fa. As shown 
by the form of the solutions, viscosity has a second-order 
effect on fa and fa, whereas it has a first-order effect on 
fa and fa. The inclusion of fa in the general solution of the 
disturbance for Blasius flow therefore indicates a greater 
effect of viscosity on the disturbance than for the case of the 
free boundary layer. As could be inferred from the fore- 
going, therefore, the effect of viscosity on the stability 
characteristics of the free boundary layer is apparent only 
at very low Reynolds number, whereas the Blasius flow 
stability characteristics are much more affected even at 
higher Reynolds numbers. 

The inaccuracies due to small values of the parameter «i? 
can be avoided by direct integration of the Orr-Sommerfeld 
equation for those cases. The asymptotic solution of the 
entire disturbance function could easily be developed as was 
done for the expanded disturbance function, and the numeri- 
cal technique of integration after starting from the asymp- 
totic solution would correspond to equation (21). These fur- 


ther solutions should be performed as soon as more high- 
speed computing-machine service can be obtained. 

CONCLUSIONS 

It is concluded that the laminar boundary layer between 
parallel streams is an unstable-flow configuration except at 
low Reynolds numbers. The method of calculation of 
stability characteristics is successful for small absolute values 
of the parameter (—i/aR). 

In comparison with Blasius type flow against a flat plale, 
instability occurs at much lower Reynolds numbers for the 
free boundary layer than for the boundary layer against a 
flat plate with no pressure gradient. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, March 21, 1949. 


APPENDIX 

COEFFICIENTS FOR ASYMPTOTIC SOLUTIONS OF AND 


a= 

1.238493 re 

«0 = 

— 1.23849316 

a i = 

1.23849316 

a«— 

-0.30962329 


0.08600647 

b l = 

a 3 
8c 


h n — 

a 3 . 

( n ) a n ' 

02 — 

c 

V 2+ 16c 1 

A — 


(—a +— 

— 

' c 1 

k 8 3+ 8c 


_1 

c 





rfl- 


i a i 


2 C 2 


hia — 


h 


■~-+G7 




_ bih w + bj 
2ao+o a 


32c 2 


ffli 


') 


20 


i ^ 1A20+ bah-iQ-\~bi 

"30 sn s J 

3 aa + (J aj 
Ki — a 2 Aio 

K 2 =(2a) ! (a+isji 

# 3 =(3a) 2 (W| aj I h0 

li=d 0 Ki 

l 2 =d 0 K i +d l K 1 

l s = d qK s -f- d i K 2 -f" d %K j 


h n 


ku — 


bihn~\~ b 2 


2aa + a 2 

l b \h%i -f- b 2 hn -{- 6 3 -j- 1 2 

fal /O \2 _ 

3aa+^2 a ) 
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TABLE I— EIGENVALUES OF a, c, AND R 


[Inviscid solution C E— »): c.=0.5S727I98; and <x,-0.39537.] 


a 

c 

R 

a. 

c 

R 

Im c, 0 

0-390 

.390 

.390 

.390 

.390 

.390 

.390 

.385 

.385 

.385 

.385 

-3S0 

.380 

.380 

.375 

.375 

.375 

.375 

.375 

.365 

.365 

.365 

.365 

.355 

.355 

0. 5S4COO 
.584550 
.584581 
.584600 
.585000 
.588000 
.600000 
.581500 
.582190 
.582204 
.582500 
.579960 
.580000 
.583000 
.575000 
.57//00 
.577850 
.578000 
.580000 
.573600 
.573900 
.574050 
.574100 
.567300 
.570000 

292.54+ 54.12i 
305. S9+ 1.961 
307.06- .011 

307.27- 1.25i 
312.08 - 28.00i 
210. 66-190.29i 
3.96- 88.131 
156.95+ 11. 59i 
160.64+ .25! 

160.76+ .011 

161.85- 5.171 
109.78- .01i 

109.83- .32i 

109. IS- 25. 381 
79.73+ 11.711 
83. 69+ . 681 

83.81+ ,01i 

83.94- .671 

84.40- 10.05i 
57.40+ ,98£ 

57.19+ .3Si 
57.52+ .071 

57.55+ .01i 

43.68+ 4.08? 
44.18+ L13i 

0.355 

.345 

.335 

.335 

.335 

.335 

.335 

.325 

.325 

.315 

.315 

.315 

.315 

.315 

.305 

-305 

.305 

.295 

.295 

.295 

.295 

.285 

.285 

.285 

.285 

3 ^ 
B . 

44. 25—0. Oli 
36.19+ .04i 
30.83+ .52i 
30.76+ . 12i 
30.74+ .021 

30.68- .371 

30. 45-1. 351 
26.92+ . 6S1 
26.75- .021 
24.35+4.661 
24.40+3.391 
24. 09+1. 291 
23.63+0 1 

22.68— 1.691 

21.12+ .161 
21.00+0 i 
21.00- .021 
18.67+ .021 
18.67+0 i 
18.63- .021 
18.56- .091 
16.38+ .101 
16.26+ .021 
18 25+0 1 

16.21- .021 

Im c, 0.05 

0.315 

.315 

.285 

.285 

.285 

.285 

0.562500+0.051 
.562800+ .05t 
.540000+ .05£ 
.555000+ .052* 
.570000+ .051 
.580000+ .05/ 

72.57+ 0.95* 
72.68- .07* 
31.19+ 8.00t 
32.55+ .2 5; 
28. 87— 7. 541 
24. 47-10.722* 

0.265 

.265 

.265 

.245 

.245 

0. 548000+0.05* 
.550000+ .05* 
.560000+ .05* 
.560000+ .051 
.570000+ .05* 

28.89+ 2.401 
24.79+ 1.901 
23.74- .791 
19.47+ 1.111 
18.26- .151 

Im c 0.10 

0. 265 
.265 
.265 
.245 
.245 
.245 

0.550000-H!.10i 
. 551700+ ,10i 
.552000+ . lOi 
.545000+ -10i 
.548000+ . IOi 
.550000+ -10i 

79.17+ 7. 14i 
81.12+ . 12i 
58.03- 1.29i 
39.83+ .881 
39.72- I.8Si 
39.45- 3.73t 

0.225 

.225 

.225 

.205 

.205 

.205 

0.542000-rQ.lOi 
.544600+ .10i 
.546000+ -10f 
.550000+ - IQi 
.551000+ .101 
.557000+ .10i 

27.74+ 1.481 
27.58+ .751 
27.32+ .031 
2L 06+ 1.091 
20.58+ .451 
20. 18+ .011 

Im c, —0.05 

0.400 

.400 

.400 

.400 

0.500000-0.052 
.550000- -052 
.580000- .052 
.590000- .05i 

16.25+14.25* 
26.91+13.76* 
34.34+ 1.64x 
33.20- 4.56* 

0.400 

.420 

.420 

0. 600000-a 05* 
.580000- .05* 
.590000- .05* 

29.59- 9.721 
48.72+10.091 
49.92- 3.101 


TABLE II— EIGENVALUES OF a, c FOR R REAL 


cc 

c 

j R 

CL 

i c 

R 

Imc, 0 

1 Im e t 0.05 

0.395 

0.587 

a 

a 3i5 

0.563+0.05* 

73.6 

.390 

.585 

307.1 

.285 

.555+ .051 

; 32.6 

.385 

.582 

160.8 

.265 

.557+ .051 

24.1 

.380 

.580 

109.8 

.245 

.569+ . 051 

18.4 

.375 

.578 

83.8 




.365 

.574 

57. 6 




.355 

.571 

44.3 


Im c, 0.10 


.345 

.569 

36.2 




.335 

.567 

30.7 




.325 

.567 

26.8 

0.265 

0.552+0.101 

81.1 

.315 

.568 

23.6 

.245 

.546+ .101 

39.8 

.305 

.571 

21.0 

.225 

.546+ .10/ 

27.3 

.295 

.576 

18.7 

.205 

.557+ .10* 

20.2 

.285 

. 5S6 

16.3 







Im c, -0.05 




0.400 

0.583 -0.05* 

34.0 




.420 

.588- .05* 

45. 9 




